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Dynamical aspects of quantum entanglement for coupled mapping systems 



m 
o 
o 

(N 



> 
in 

o 

O 
m 
o 

Ok 
^— > 



X 



Hiroshi Fujisaki*, Atushi Tanaka 1,1 ', and Takayuki Miyadera 2 ^ 

Department of Theoretical Studies, Institute for Molecular Science, Myodaiji, Okazaki, 444-8585, Japan 
1 Department of Physics, Tokyo Metropolitan University, Minami-Osawa, Hachioji 192-0397, Japan 
2 Department of Information Sciences, Tokyo University of Science, Noda City, 278-8510, Japan 

(Received February 1, 2008) 

We investigate how the dynamical production of quantum entanglement for weakly coupled 
mapping systems is influenced by the chaotic dynamics of the corresponding classical system. 
We derive a general perturbative formula for the entanglement production rate which is defined 
by using the linear entropy of the subsystem. This formula predicts that the increment of the 
strength of chaos does not enhance the production rate of entanglement when the coupling is 
weak enough and the subsystems are strongly chaotic. The prediction is confirmed by numerical 
experiments for coupled kicked tops and rotors. We also discuss the entanglement production 
using the Husimi representation of the reduced density matrix. 
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1. Introduction 

Quantum information processings 1 '' have been recog- 
nized as a new paradigm of science, because of its fun- 
damentally, applicability, and interdisciplinarity. For a 
quantum computer, a typical example of the paradigm, 
to be really effective, it must consist of a large number of 
(interacting) qubits, and its size cannot be microscopic. 
Such a "complex" structure of an effective quantum com- 
puter causes many problems due to decoherence 2 ) and 
"quantum chaos" . 3 - ) 

Here we investigate the relation between (quantum) 
entanglement (which is an essential ingredient of quan- 
tum information processings) and chaotic behavior of the 
corresponding classical systems. The systems we employ 
are weakly coupled mapping systems, and each subsys- 
tem can be chaotic in the classical limit. We can imag- 
ine as follows: A quantum computer has many process- 
ing units which are "chaotic" in some sense (e.g., in the 
classical limit), and they weakly interact with each other. 
How much entanglement is generated between the units 
in such a situation? This question is interesting from the 
view point of robustness of quantum computation. 3 -* 

Due to the weakness of the coupling, we utilize a per- 
turbation theory to analyze the entanglement production 
in the system considered. In this paper, we discuss gen- 
erality of our formula by employing two typical chaotic 
systems, i.e., coupled kicked tops and rotors (Sec. 2). Nu- 
merical results for both systems are presented (Sec. 3). 
To complement our previous papers, 4 ' 5 - ) we also discuss 
the wavefunction properties of the subsystems in the en- 
tanglement production region employing the Husimi rep- 
resentation (Sec. 4). 

2. Coupled mapping systems and the perturba- 
tive formula for entanglement production 

Here we consider a composite system which consists 
of two subsystems. We denote one-step time-evolution 
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operators for each subsystem as U\ and U%. We also 
introduce another one-step time-evolution operator U e 
which describes the interaction between them. Further- 
more, the coupling time-evolution operator is assumed 
to be 



U e = exp{-ieV/h} 



(1) 

where V = 8> Qa ^ with q$ being a a-th dynam- 

ical variable for subsystem i, and e is a coupling param- 
eter. The latter is used in the perturbative treatment 
below. 

Hence the one-step time evolution for the whole system 
is described by 



\*(t+l)) = U e U 1 U 2 Mt)). 



(2) 



Since we only examine the case where the whole system is 
in a pure state, we quantify the entanglement production 
by the linear entropy of the subsystem 



SUt) = i-Tn{p {1) (t) 2 }, 



(3) 



where p^\t) is the reduced density operator for the first 
subsystem. 

Using the time-dependent perturbation theory, we de- 
rive a formula (see Refs. 4, 5 for a simpler case and its 
derivation) for the linear entropy: Su n (t) = S^(t) + 
C(e 3 ) as 



t t 



(4) 



S^(t) = S ^^D(l,m) 

1 = 1 m=l 

where So — 2e 2 /h 2 and D(l,m) is a correlation function 
of the uncoupled system: 



D(/,m)=^C'W(/,m)C'3(/, i 

a, 13 



(5) 



with 

C«, (hrn) = (^(O^H), - (?W(0>i<$V)>< (6) 

and qa (I) = {U[)^qit U[ etc. represents a free time evo- 
lution of the subsystem's variable without interaction. 
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In the next section, we apply this result to two model 
mapping systems, i.e., coupled kicked tops and rotors. 

3. Stronger chaos does not imply larger entan- 
glement production rate 

We introduce two coupled kicked systems as follows. 
Coupled kicked tops are described by the following uni- 
tary operators: 



Uy = 
U 2 = 

u f = 



e -ikiJlj{2jh) e -iKj yi /(2ft) 



e -ik 2 Jl 2 /(2jh) e -iTT.J V2 /(2ft) 



-ieJ zl J z2 /(jfi) 



and coupled kicked rotors are by the following: 
U 2 



p-ilf I \2K) g-ifci cos^i/ft 
p -ili I (2ft) p -ik 2 cose 2 /fi 



u e = 



-2ecos(#i— &2)/h 



(7) 

(8) 
(9) 

(10) 

(11) 
(12) 



where fci, k 2 represent the strength of nonlinearity re- 
lated to chaotic properties of the systems. For the de- 
tailed information of the former and latter systems, see 
Refs. 4, 5 and Ref. 6, respectively. We note that the pe- 
riodic boundary conditions with period 2ir are imposed 
for both variables 0, and /j. 

As an initial state, we take a product, i.e., separable 
state of spin coherent states 5 - 1 (coherent states) for cou- 
pled kicked tops (rotors). In Fig. 1, according to Eq. (3), 
we calculate the time evolution of the linear entropy for 
the "chaotic" initial conditions, i.e., their corresponding 
classical states in phase space are embedded in chaotic 
seas. For both cases, the linear entropy increases lin- 
early as a function of time in this transient region, which 
is called i-linear entanglement production region in the 
following, before the "equilibration" of the entropy. In 
these cases, the coupling is in a sense weak, and we apply 
the formula, Eq. (4), to these situations. 

From a classical intuition, we expect that the correla- 
tion functions D(l,m) decay very quickly as a function 
of \l —m for strongly chaotic cases (this is numerically 
confirmed for a kicked top 4,5 ) and a kicked rotor 7 )), so 
it is assumed to be 



D(l,m) ~ D Q e-^ l - m \. 



(13) 



Using above, we can easily derive the following for the 
entanglement production rate: 



jcPT 
"^lin 



(*) 



dt 



t»l/7 



r coth( 7 /2) 



(14) 



with r = SqD . Since 7 often becomes large as the 
nonlinear parameters increase, and there is a numerical 
experiment that 7 and the sum of positive Lyapunov 
exponents are correlated in the coupled kicked tops, 5 ) 
this formula implies that strong chaos (7 — > 00) leads to 
a saturation of the entanglement production rate (r — ► 
r ). This prediction is actually confirmed by numerical 
experiments for both systems (coupled kicked tops and 
rotors) as shown in Fig. 2. 
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Fig. 1. Time evolution of the linear entropy for the coupled kicked 
tops (a) and rotors (b). The initial states are embedded in 
chaotic seas in both systems. For simplicity, we take k\ = 
&2 = k. The interaction strength is weak (e = 10~ 4 for tops 
and e = 10 -3 for rotors), ft = 1.0 and j = 80 for tops and 
ft= 2tt/128 for rotors. 



The above is a result for strongly chaotic cases. In 
contrast, we comment for weakly chaotic cases using the 
above formula: For strongly chaotic cases with bounded 
phase space, Do saturates to a certain value determined 
by a statistical argument, whereas, for weakly chaotic 
cases, D grows up as the nonlinear parameters increase. 
This roughly explains the common belief: chaos en- 
hances entanglement for weakly chaotic systems. (Sec 
the cited references in Refs. 4,5.) This seems to contra- 
dict our result, but the point is that we focus on strongly 
chaotic cases where D saturates, and in such a situa- 
tion, stronger chaos saturates the entanglement produc- 
tion rate for weakly coupled mapping systems. We also 
comment that this formula can be easily extended to 
describe flow systems with continuous time, 5 ) which is 
useful when applying to more realistic systems. 

4. Husimi representation as a tool to study en- 
tanglement production 

The linear entropy which we employ as a measure of 
entanglement is calculated from the reduced density op- 
erator p( x \t). Hence it is natural to ask what happens 
in p^(t) itself when the entanglement production oc- 
curs. Here we shall address this issue using the results 
of coupled kicked tops. 
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Fig. 2. Nonlinear parameter dependence of the normalized entan- 
glement production rate (T/To) for tops (a) and rotors (b). The 
plots for five different initial conditions embedded in the chaotic 
seas are superimposed. We take k\ = ki = k. The k depen- 
dency of r for the coupled rotors exhibits oscillations. We sug- 
gest that the oscillations are similar ones that arc observed in 
the diffusion constant of a kicked rotor due to higher-order time 
correlations. 8 ) 




Fig. 3. The contour plot of the absolute value of (a) the density 
matrix at t = 112 for a kicked top, and (b) the reduced den- 
sity matrix t = 112 for coupled kicked tops with e = 0.003. 
The initial state is embedded in the chaotic sea, and given by 
\0,<j>) = |0.89,0.63) for (a) and \0 U 4>i) ® \8 2 , <h) = |0.89, 0.63) ® 
|0.89,0.63) for (b). 



In Fig. 3 (b), we show the absolute value of the reduced 

density matrix Pm\,m 2 (t) — 0' TO i|P^(*)U TO 2) during the 
^-linear entanglement production region (see Sec. 3). We 
also show the corresponding density matrix of a single 
kicked top [Fig. 3 (a)]. Here \jm) is the simultaneous 
eigenstate of J 2 and J z for top 1. From the previous 
works, 11 ) we expect that the decay of the off-diagonal 
elements of the reduced density matrix reflects entangle- 
ment production. Such an analysis is useful when the 
entanglement production becomes almost equilibrated, 
however, two cases in Fig. 3 are very hard to distinguish 
since the system is in the t-lincar entanglement produc- 
tion region. (Of course, after a long time, the off-diagonal 
elements finally decay.) To solve this issue, we propose 
to examine the Husimi representation of the density op- 
erator. 12 ) 

The Husimi function for the kicked top system is de- 
fined by 

iT(M) = <M|p (1) (*)IM>, (15) 

where \6, <p) is a spin-coherent state given by 

(jm\0, <t>) = T-, i ioV- \ / T- 777 1 77 (16) 

(1 + |7| 2 ) J V {J +m)\(j -m)\ 



with 7 = e 1 ^ tan(0/2). For simplicity of visualization, we 
take a smaller value of j (j = 30 instead of j = 80 used 
in Figs. 1 and 2). Though the normal-scale plots do not 
show any significant fingerprint of entanglement produc- 
tion due to the interaction (Fig. 4), the log-scale plots do 
show it (Fig. 5). The Husimi zeros are easily confirmed 
in Fig. 5 (a). They become local minima with positive 
values as the interaction strength gets larger [Fig. 5 (b)]. 
Husimi zeros are considered to characterize "complex- 
ity" of the wavefunction, 9 ) so the quantification of this 
local minima, the traces of the Husimi zeros, might char- 
acterize the entanglement production from "complexity 
theory" point of view. 10 ) Many measures for complex- 
ity of a Husimi function have been proposed so far, 9, 12 ' 
but we suppose that another measure, e.g., the phase- 
space averaged curvature of a Husimi function, is needed 
to quantify the complexity related to entanglement pro- 
duction. This issue will be discussed elsewhere. 

5. Summary and outlook 

We have investigated how the production rate of quan- 
tum entanglement is affected by the chaotic properties of 
the corresponding classical system using coupled kicked 
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Fig. 4. The contour plots of the Husimi functions corresponding 
to the density matrices in Fig. 3. The scale is normal. 




Fig. 5. The contour plots of the Husimi functions corresponding 
to the density matrices in Fig. 3. The scale is logarithmic. Note 
that the number of the Husimi zeros is 2j = 60 [see (a)]. 

tops and rotors. We have derived and numerically con- 
firmed that the increment of the strength of chaos does 
not enhance the production rate of entanglement when 



the coupling is weak enough and the subsystems are 
strongly chaotic. We believe that this conclusion is gen- 
eral and applicable to any mapping system with bounded 
phase space. It is often believed that chaos enhances en- 
tanglement, but that can be the case for weakly chaotic 
regions, and such effects can be also explained by the 
factor Do in our formula, Eq. (14). Note again that the 
above our conclusion is basically for strongly chaotic re- 
gions where D saturates. We also proposed to use the 
logarithmic plot of the Husimi representation of a re- 
duced density operator to investigate entanglement pro- 
duction in coupled quantum systems. 

Since the dynamical aspects of entanglement produc- 
tion have been studied in relation with quantum com- 
puting, 13 ) we hope that our analysis will be useful in the 
studies of quantum information processsings. In Rcf. 3, 
Prosen and Znidaric showed that quantum computing 
can be more robust with the use of quantum chaotic 
systems compared to the corresponding regular systems. 
We expect that the saturation of entanglement produc- 
tion found here will be utilized in such a situation. In ad- 
dition, entanglement production is strongly related to de- 
coherence processes, 2 ^ so we expect that there are many 
applications of our results for quantum dynamical pro- 
cesses like quantum optical processes in environments 14 - ) 
or chemical reactions in solvents. 15 ) 

One of the authors (H.F.) thanks T. Prosen, K. Taka- 
hashi, M. Toda, S. Kawabata, and K. Saito for useful 
comments. 
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